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Abstract 

We consider here the bremsstrahlung emission of photons at low and intermediate energies 
Eiab — 1000MeV/u of the projectile, and derive expressions more general than previous 
results obtained by Neuhauser which were limited to the case of isotropic systems. We 
find that the two-photon correlation function strongly depends not only on the space-time 
properties of the collision region but also on the dynamics of the proton-neutron scattering 
process in nuclear matter. As a consequence of polarisation correlations it turns out that 
for a purely chaotic source the intercept of the correlation function of photons can reach the 
value 3 (as compared with the maximum value 2 for isotropic systems). Furthermore even 
for "hard" photons [E 1 > 25MeV) the maximum of the correlation function can reach the 
value of 2 in contrast with the value of 1.5 derived by Neuhauser for this case. The formulae 
obtained in this paper which include also the possible presence of a coherent component can 
be used as a basis for a systematic analysis of photon intensity-interferometry experiments. 



* E . Mail : razumov@convex.hrz . Uni-Marburg. de 
f E. Mail: WEINER@vax.hrz.Uni-Marburg.de 



1 Introduction 



The Hanbury-Brown and Twiss intensity interferometry Q plays at present an important 
role in particle and nuclear physics in a wide range of cm. energies, because this method 
can provide information about the space-time properties of the sources. From hadron inter- 
ferometry one obtains information mainly about the properties of the fireball just near the 
freeze-out. Moreover, in this case we are faced with some dynamical problems like the final 
state interaction and contribution from resonances which mask the real space-time size of the 
collision region. In contrast, photons and leptons have a large mean free path and leave the 
collision region just after their production. Additionally we can assume that photons have 
no final state interaction (light-light scattering can be neglected). Therefore direct photons 
(i.e. photons not originating from 7r° decay) and dileptons are good probes for the earliest 
stages of the reaction. Another advantage of electro-magnetic probes is that their emission 
is controlled by QED rather than by strong interaction dynamics for which no comparable 
theory exists yet. 

In this paper we consider mostly the emission of "hard" photons (their definition will be 
given below) in heavy ion reactions in the energy range up to 1 Gev/u. We will show that 
the photon-photon correlation function can provide there not only space-time information 
but also some dynamical information about the proton-neutron scattering in nuclear matter. 
It is known |2], || that at intermediate energies hard photons are produced mainly through 
bremsstrahlung and 7r decay. Since the photon pairs from the ir Q decay can be eliminated 
to a great extent by measuring the photon pair invariant mass, we concentrate ourselves on 
bremsstrahlung photons. 

2 Bremsstrahlung photons and their correlations 

The photon spectrum is usually described by a superposition of two contributions, both 
parametrized by exponentials [|2], |3| and referring to different photon energies respectively: 
the low energy part is assumed to be due to thermal photons and the slope is just given by 
the temperature. The second (high energy) part can be considered || as due to a mixture 
of bremsstrahlung (see Chapter |^) and photons originating from decays (mostly ttq — > 77). 
Measurements of the angular distribution of radiation, the photon-source velocity 0, |3[] as 
well as the impact-parameter dependence of the photon multiplicity || suggest that the 
direct hard photons mainly originate from bremsstrahlung in independent proton-neutron 
collisions^. 

As was shown in || the non-relativistic recoilless bremsstrahlung formula for the current 

lr The photon emission from proton-proton collisions is of quadrupole form and therefore highly suppressed 
as compared with the dipole radiation from proton- neutron collisions. 
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(where p = pj — pf is the difference between the initial and final momentum of the proton, 
e\(k) is the vector of linear polarisation and k the 4-momentum) works well even in the 
relativistic case due to the fact that relativistic and recoil corrections to some extent com- 
pensate each other. If we have N such sources we can write down the transition current as 
follows: 
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The index n here labels the independent p — n collisions taking place at different space- 
time points y n . Formulas flU), (0) are examples of radiation by classical currents, where the 
influence of the emitted photons on the p — n collision process is negligible. 
For an arbitrary classical current J^{x) coupling with the photon field 
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through the conventional interaction lagrangian 
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we can easily find the exact S'-matrix 
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which contains all the information on the photon production and absorption. We use in ([5]) 
the notation: 

j\k) = ^(k)r(k) . (6) 

where J^(k) is the Fourier transform of the current J^(x). Before the collision there are no 
photons, i.e. we have the photon vacuum in the initial state. 

The amplitude to produce n-photons from the photon vacuum by the classical current J 
follows directly from (pT): 
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2 The experimently found one-photon spectra are of exponential form This property can be also 

incorporated into this model, see Chapter ^. 



The corresponding cross-section can be written in a condensed and convenient form by 
means of the generating functional for the exclusive processes. Let us introduce the notation 
q = (k, A) for the momentum k and polarisation A degrees of freedom of photon. Integration 
over q means integration over momentum (with the invariant measure dk ) and summation 
over polarisation A . We define the exclusive generating functional g ex [Z] as follows: 
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The conservation of probability gives the normalisation g ex [Z = 1] = 1 . The method of 
the generating functionals allows one to avoid the complicated combinatorics which usually 
appears when one consider inclusive processes. The generating functional for the inclusive 
cross-section can be introduced in a similar way 
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(with the normalisation condition g m [Z = 0] = 1 ) and can be very simple expressed through 
the exclusive one 

g in [Z] = g ex [Z + l] . (10) 

In our case (see eq.(^)) the summation over n can be performed explicitly and we are left 
with the simple expressions 
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If the classical current obeys some random behaviour which is the present case, as we consider 
chaotic sources, the generating functional should be subjected to the averaging over the 
current distribution 

G in = (g*n\Z^ . (13) 

Taking the variational derivatives of G m [Z] with respect to Z one gets the inclusive spectra. 
For example single and double inclusive spectra read as follows: 
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An essential ingredient in eqs. ( pT3|) - (|T5| ) is the average prescription < ... >. 

For the model considered here, where hard bremsstrahlung photons originate mostly 
from independent proton-neutron collisions at the early stage of the heavy-ion reaction, the 
averaging prescription < ... > appears quite naturally and consists of two parts: 



1. The points y n where the independent p — n collisions take place (see eq. (0)) are 
considered to be randomly distributed in the space-time volume of the source (e.g. 
fireball) with a distribution function f(y) for each p — n collision. This kind of 
averaging is typical for Bose-Einstein correlation studies. 

2. The amplitude of photon production ([!]) is very sensitive to the momentum transfer of 
the proton in the proton-neutron collision. That is why we have to take into account 
also the fluctuations of the initial momentum of the proton: p» = po + Ap^ . Here po 
is the momentum of the nucleus as a whole and Ap^ is the Fermi- motion of a nucleon. 
The final momentum distribution of the proton (pj) can be in principle determined 
by dynamical models of nucleus-nucleus collisions, e.g. BUU ||, etc. 

Since in practice we need only the first two moments of this distribution we keep them as a 
free phenomeno logical parameters to be determined in the experiment. The comparison of 
the experimental values with model predictions can serve as an indirect test of the model. 

The fundamental quantity in our approach is the two-current correlator in momentum 
space: 
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- f(h)f(-k 2 ) < fcfa) >< j x *(-k 2 ) >} + < J Al (^i) >< J A2 (-A: 2 ) > . (16) 

Here we use the properties (J x (k)J = J x (-k), < j Xl (k 1 )j^(-k 2 ) >=<j Xl (h) >< j%{-k 2 ) > 
for n m (that is the hypothesis of independent proton-neutron collisions), and f(k) de- 
notes the Fourier transform of f(y) with the normalisation f(k = 0) = 1 . 
Usually the function f(k) has a steep maximum around k = with a width of the order of 
the size of the source R . In the region of the Bose-Einstein peak for hard photons in partic- 
ular for 2kR > 1 we can therefore neglect f(ki)f(—k 2 ) as compared to f{k\ — k 2 ) . This 
approximation strongly simplifies the algebra without influencing significantly the accuracy 
of the calculations. 

3 Chaotic sources 

Let us consider the case of a chaotic source for which < J X h{k) >= . The current correlator 
(|16D is now given by the simple formula: 
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< p^p^ > denotes here the averaging with respect to the distribution of the momentum 
transfer of the proton in the collision n and Y,n=i goes through all the relevant proton- 
neutron collisions. As mentioned above the quantity < p^p{ > can be extracted from 
dynamical models of heavy-ion collisions. But one can reach important and general conclu- 
sions without specifying this quantity as follows. We use the axial symmetry around the 
beam direction. The tensor decomposition of < p^p{ > gives then: 

< PnPi > = + ^ ' (18) 

where £ is the unit vector in the beam direction and a n , 5 n are real positive constants. Note 
that this expression is more general than the corresponding one used in || where because 
of the isotropy assumption 5 n was assumed to vanish. This generalization has important 
consequences to be exhibited below. Now let us separate the average value of the proton 
momentum transfer < p n > : 

p n = Ap n + < p„ > (19) 

We have then: 

< PnPi >=< ^PlAK > + < pj, >< pi >=< Ap^Ap* > + < p > 2 €V (20) 

where the tensor < Ap^Ap^ > can be represented (due to axial symmetry) again as the 
sum of the two terms: 

< Ap^ n > = ^ + t^et? . (21) 

In order to find the coefficients a n , £ n , 5 n we split the 3- vectors in the transverse and the 
longitudinal parts: 

Ap l n = £-(£-A Pn ); Ap* n = A Pn - Ap^ . (22) 



The coefficients in fljjf), fl2T|) read 
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With the help of (|l8|) , fl23"l) one can express the current correlator (|l7|) 
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through two parameters a and 5 



a 



i " i 



t \2 



> 



1 



N 



iV n=l 



a 

3 ' 



(25) 



(26) 



which absorb the relevant dynamical information about the proton-neutron scattering in the 
medium. The parameters a and 5 can be extracted from experimental data on photon- 
photon Bose-Einstein Correlations (cf. below) and/or calculated from dynamical models for 
heavy-ion collisions. 

The single inclusive cross-section for a detector, which does not measure polarizations 
follows directly from ([14]), ([24]) 
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where 9 is the angle between the photon and the beam directions and the polarization sum 
is calculated using the well-known identity: 
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with n = k/|k| . 

To calculate the double inclusive cross-section (as well as the higher order inclusive spec- 
tra) one has to know higher order current correlators. In our case when all proton-neutron 
collisions are considered to be independent from each other and the number of the partic- 
ipating protons is sufficiently large N > 10 , we can apply the central limit theorem and 
express the higher order current correlators through the first and second ones. Assuming 
< J x {k) >= (no coherence) the double inclusive cross-section is represented through the 
sum of products of the two-current correlators 

P2(h,k 2 ) = Yl <J Xl (h)J Xl (-k 1 )J x i(k 2 )J x *(-k 2 )>= 
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and the polarization sum 
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is performed with the help of (|24]), (|28|). 

The general expression for the second order correlation function is defined by 

C 2 {ki,k 2 ) = 31 
pi(fci)pi(fc 2 ) 

and has a complicated angular dependence, but the two limiting cases 1) a ^> 5 and 2) 
a C 5 lead to very simple expressions: 

1. For the case a > 5 we have: 

C 2 (h,k 2 \a^ 0,6 = 0) = 1 + ^(1 + cos 2 ^)[|/>i-A;2)| 2 +|/(A;i + A;2)| 2 ] (32) 
which is the result derived in || and which gives for the intercept 

C 2 (k,k) = 3 -+ 1 -\f(2k)\ 2 . (33) 

2. The opposite case a <C 5 leads to another formula: 

C 2 (k h k 2 \o = 0, 6^ 0) = 1 + \f(h - k 2 )\ 2 + \f(ki + k 2 )\ 2 ■ (34) 
with the intercept exceeding 2: 

C 2 (k,k)=2 + \f(2k)\ 2 . (35) 

For hard photons \f(2k)\ 2 1 (cf. ||) and one can neglect this contribution to the 
intercept while for soft photons f(2k) is non-neglegible and in the limit k = 0, f(0) = 1 so 
that C 2 (k,k) = 3. The real situation (when both a and 5 contribute) is between ( j32f ) and 
( |34D and exhibits a more complicated angular behaviour than the considered above limiting 
cases 0. For instance the general expression for the intercept which follows directly from 
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depends on both parameters a and 5 and varies in the range between 3/2 and 3. The solid 
angle integration over all possible orientations of the photon momentum k corresponds to a 
47r detector. 

The function f(x) and its Fourier transform f(k)) reflects the space-time properties of the 
photon source. It depends in principle on three constants: the time duration R , the longi- 
tudinal radius Ri and the transverse radius Rt. One can propose a concrete parametrization 



for the source geometry f(k) and then fit the inclusive data using ([H]). However, because 
of insufficient statistics one has usually to limit oneself to a smaller number of parameters. 



In particular the choice of two parameters T = Rq and R = Ri = Rt is good enough for 
the present state of the art. With this assumption we can perform analytically the angular 
integration over 9\ and 9 2 in (p7|) and (p9|), (p0|) keeping constant the angle between the two 
photons rix • n 2 = cosip = const. The corresponding integration 

J dji = 2tt J sinQ\dQ\ J dip (37) 


extends over the rotations of n 2 around ni {dip) and over all orientations of ni around the 
beam direction. After some algebra we are left with the general expression for the Bose- 
Einstein correlation function 
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which has a pronounced dependence not only on the space-time characteristics Rq, R but 
also on the angle ip and the "dynamical" constants a, 5. The limiting cases (^) and ( p4|) 
as well the intercept formula (|3T)|) can be rederived directly from (133) . 



4 Partially coherent sources 

In general, photon sources are not totally chaotic but may contain also a coherent compo- 
nent < J x {k) >= I x (k) 7^ so that the total current = J x h + I x leads to a partially 
coherent field. There are several mechanisms which are responsible for coherence: collective 
deacceleration of the initial nuclei, collective flow, coherent radiation from nuclear fragments, 
etc. The possibility to investigate such collective phenomena is an interesting and important 
part of heavy-ion physics. In this chapter we discuss phenomenologically the influence of 
the coherent part of the electric current on the photon correlation function. The single and 
double inclusive cross-sections (see (|27|) and (^) ) read: 

Pi(k) = j:(F x \k,k) + \I x (k)\ 2 ) , (39) 
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As before we assume axial symmetry around the beam direction and parametrise I x (k) as 
follows: 
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Using (H, (|D and fl41]) we get 
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(43) 



which together with F xx (k,k) and F XlX2 (ki, k 2 )F XiXl (k 2 , k x ) (see (|27J) , (|0]) ) determine com- 
pletely the single and double inclusive spectra (^,0). As compared to the completely 
chaotic case there is now one more function S(ko) entering in pi(k) and p 2 {ki, k 2 ). 

Significant simplifications can be achieved again assuming Ri = R t = R. The time 
dependence of the photon source is still arbitrary. It is then possible again to perform 
explicitly the integrations over all the angles except for ip and write the photon-photon 
correlation function in the compact form 

c 2 (h, k2 ) s /,^f;- fc ;> = (44) 
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= 1 + {\ 2 (k )A\f(k 1 -k 2 )\ 2 + 2\(k )(l-\(k ))B'JZef(k 1 -k 2 ) + (k 2 ^-k 2 )} 
where we have introduced a new set of parameters: 

a + b + \S(k )\ 2 a + b 

1 1 + cos 2 ip + (1 -xf (13 + cos 2 ip)/ '5 
~4 1 + (1 - x\) 2 (2>cos 2 il) - 1)/10 ' 

1 1 + cos 2 ip + (1 - x) (13 + cos 2 ^) /5 
~4 ' 1 + (1 -xA) 2 (3cos 2 V» - 1)/10 ' 
Here A(fco) is of the chaoticity parameter which is the ratio of the number of chaotically 
produced photons with energy k and their total number A(fco) =< n(k ) > c / l / < n(k ) > . 
We assume also that on the scale important for the Bose-Einstein study (|fei — k 2 \ ~ 1/-R) 
the function X(k ) does not vary significantly: - k%)d\(k°)/dk°\ < (A(fc°) + A(fc§))/2 . 
The other parameter which influences strongly the behaviour of the correlation function 
C 2 (ki,k 2 ) is the isotropy parameter x (see ([45])). Both these parameters A and x can be 
extracted from the analysis of the angular distribution of the radiation and the value of the 
intercept C 2 (k, k): 

Pl (k)/ <n(k )> = ^L[ x \(k ) + l(l-x\(k ))sin 2 9] , (46) 
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C2 ^ k ) = 1 + ^+{\ 2 -l\)* { A2(A:0)[5 + 7(1 _ ^ + 2A(1 " A)[5 + 7(1 " X)] ) (47) 



(for the case of hard photons |/(2A;)| 2 <C 1 and one can neglect this contribution in (|47|) ). 
After this has been done the only unknown function in (|44"D is the space-time distribution 
function of the source f(x) (or f(k)) which can now be obtained by fitting the experimental 
data. This function has the physical meaning of the space-time distribution of the radiating 
region and reflects the geometry of the early stages of the collision. 

We would like to stress once again that both the parameters A and x strongly influence 
C^&i, k 2 ) and the knowledge of only the one of them (e.g. A) is not enough to determine in 
a unique way the two-photon correlation function. For instance, one can check that all the 
A, x pairs, for which A = (y/2 - l)/(y/2 - x) and < x < 1 , lead to the same intercept 
C 2 (k,k) =3/2 . 



5 Exponential fall of the bremsstrahlung amplitude 

The experimental observations ^ |||] show that the one-particle inclusive spectrum of brems- 
strahlung photons has exponential form suggesting that the underlying proton current reads: 
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rather than ([[J). The single and double inclusive cross-sections calculated with (f|8|) instead 
of (P can be obtained from the previous results multiplying them by the corresponding 
powers of exp[-\k°\/{2E )} and exp[-\k2.\/(2E )} . 



The homogeneous functions like the two-photon correlation function ([H]) and the angular 
distribution of radiation fl46| ) remain unchanged and all the conclusions about the Bose- 
Einstein correlations obtained above hold. 

In the following we shall derive the non-relativistic classical trajectory of a proton which 
leads to the mentioned above current ( fi8|) . 

The current in momentum space is defined through the proton trajectory as follows: 
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(49) 



In non-relativistic case \v(t)\ 1 and as a consequence kot — kr(t) = k$t . Therefore j x (k) 
reduces to the time Fourier transform of the velocity. Using the identity 
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one finds the trajectory 
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leading to the proton current (fHf) . The standard formula (|l]) corresponds to the special case 
E — > +oo when the proton trajectory is described by the step function v(t) = v Q(—t) . 



The finite value of Eq reflects more smooth then step-like deacceleration of the proton with 
the characteristic stopping time r ~ 1/ (2E ) and stopping length I ~ vqt. For instance, for 
the projectile energy 45 MeV/u (slop-parameter Eq = 18 MeV ||) one gets r = 5.6 fm/c 
and I = 1.7 /m which seem to be quite reasonable. 



6 Summary 

In this paper the production of photons is analysed in the framework of quasi-classical 
approximation. Our consideration is valid in the region of low and intermediate colliding 
energies up to 1000 MeV ju. We assume that the hard photons (E 1 > 25MeV) are produced 
in independent proton-neutron collisions^ (see Chapter 2) and the whole system has the 
axial symmetry with respect to the beam direction. Then we derive the expressions for the 
single and double inclusive cross-sections and for the two-photon correlation function as well 
(see Chapters 3,4). We show that the behaviour of the photon-photon correlation function 
depends not only on the space-time properties of the collision region (function f{k\ — k 2 )) 
but also on dynamics of the proton-neutron scattering in matter (parameters a, 5 and the 
amount of chaoticity A). So far the photon intensity interferometry can be considered as an 
indirect way to check the dynamical properties of the heavy-ion system as well as to obtain 
the space-time information about a collision. 

It turns out that the maximum value of the correlation function (the intercept C 2 (k, k)) is 
very sensitive to the details of the proton-neutron scattering and varies generally speaking in 
the interval 1 < C 2 (k, k) < 3 (see ©). If we deal with so hard photons that \f(2k)\ 2 < 1 
the intercept finds itself in more narrow range 1 < C 2 (k, k) < 2 (which is nevertheless wider 
then the one obtained in || 1 < C 2 (k, k) < 1.5 ). The reason why the intercept can exceed 
the value 1.5 (see |J) is strongly connected with the beam-direction anisotropy specific for 
heavy-ion collisions. 

The photons produced with the close momenta can have the correlations in their polarisa- 
tions. As far as the photons with the same polarisation obey the Bose-Einstein effect like 
the scalar bosons (e.x. n°) and ones with the perpendicular polarisations behave like non- 
identical particles, these correlations in polarisations affect positivly on the Bose-Einstein 
peak increasing intercept. The presence of coherence leads to decreasing of the correlation 
effect. In order to study the two-photon correlation function including the influence of the 
anisotropy and the coherent contribution as well ( 44[) one needs additional information which 



can be obtained from the angular distribution of the radiation ( |4"6| ) and the intercept value 
fl4"7]) analysed together in order to extract the values A and x. After it has been done the 
only unknown quantity in (|44]) is the Fourier transform f(k) of the space-time probability 
distribution f(x) of the photon source which reflects the geometry of the early stages of the 

3 the photon emission from the proton-proton collisions has the quadrupole nature and therefore is highly 
suppressed with respect to dipole radiation from the proton-neutron collisions. 



heavy-ion collision. 
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